Extraction of Plumes in Turbulent Thermal Convection 
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We present a scheme to extract information about plumes, a prominent coherent structure in 
turbulent thermal convection, from simultaneous local velocity and temperature measurements. 
Using this scheme, we study the temperature dependence of the plume velocity and understand the 
results using the equations of motion. We further obtain the average local heat flux in the vertical 
direction at the cell center. Our result shows that heat is not mainly transported through the central 
region but instead through the regions near the sidewalls of the convection cell. 

PACS numbers: PACS numbers 47.27.-i 



The Rayleigh-Benard convection system consists of a 
closed cell of fluid heated from below and cooled on the 
top. The equations of motion, in Boussinesq approxima- 
tion, are 0: 

^ + v ■ Vv = - V> + iA7 2 v + gaSTi (1) 

<|E + v • VT = kV 2 T (2) 
V ■ v = (3) 

where v is the velocity field, p the pressure divided by 
density, T the temperature field, and z is the unit vec- 
tor in the vertical direction. Furthermore, ST = T — Tq 
where T is the mean temperature of the bulk fluid, g 
is the acceleration due to gravity and a, v, and K are 
respectively the volume expansion coefficient, kinematic 
viscosity and thermal diffusivity of the fluid. The state of 
fluid motion is characterized by the geometry of the cell 
and two dimcnsionlcss parameters: the Rayleigh number, 
Ra = agAL 3 /(vK,), which measures how much the fluid 
is driven and the Prandtl number, Pr = v/k, which is 
the ratio of the diffusivities of momentum and heat of 
the fluid. Here A is the maintained temperature differ- 
ence between the bottom and the top, and L is the height 
of the cell. When Ra is sufficiently large, the convective 
motion becomes turbulent. 

In turbulent convection, local velocity and tempera- 
ture measurements taken at a point within the convection 
cell display complex fluctuations in time. On the other 
hand, visualization of the flow reveals recurring coherent 
structures. One prominent coherent structure is a plume, 
which is a mushroom-like flow generated by buoyancy. 
Thus at least two strategies can be employed to study 
turbulent thermal convection or turbulent flows in gen- 
eral. One is to analyze and understand the fluctuations of 
the local measurements. The other is to characterize the 
coherent structures and study and understand their dy- 
namics. These two approaches are not independent but 
provide complementary knowledge of turbulent flows. In 



particular, there is the natural question of whether and 
how information about the coherent structures can be 
extracted from the local measurements. 

For turbulent flows not driven by buoyancy, vari- 
ous methods including proper orthogonal decomposition, 
conditional sampling and wavelet analysis have been pro- 
posed to identify coherent vortical structures from lo- 
cal velocity measurements 2j. On the other hand, much 
less work has been done in identifying plumes or ex- 
tracting information about plumes in turbulent thermal 
convection 0, Q. Belmonte and LibchaberQ used the 
skewness of the temperature derivative as a signature 
of the plumes. Zhou and Xia£| associated the differ- 
ence in the skewness of the positive and negative parts of 
the temperature difference with the presence of plumes 
and identified the plumes whenever the temperature dif- 
ference becomes larger than a chosen threshold 0. In 
Ref.[5|, plumes are identified when each of the tempera- 
ture, vertical velocity or vertical vorticity is larger than 
some threshold. 

In this Letter, we present a scheme to extract informa- 
tion about plumes using simultaneous local velocity and 
temperature measurements. Our method makes explicit 
use of the physical intuition that plumes are generated 
by buoyancy and thus the velocity of a plume should 
correlate with the temperature fluctuation. Using this 
method, we obtain the temperature dependence of the 
plume velocity and the average local heat flux in the ver- 
tical direction at the cell center. 

The experimental measurements that we use to il- 
lustrate and test our scheme of plume extraction were 
taken in an aspect-ratio-one cylindrical cell of height L — 
20.5 cm and filled with water 7]. The velocity v(t) was 
measured using a two-component laser Doppler velocime- 
try (LDV) system^ while the temperature T(t) was 
measured using a small movable thermistor of 0.2 mm 
in diameter, 15 ms in time constant, and 1 mK/f2 in 
temperature sensitivity. The velocity and temperature 
measurements were simultaneously taken using a multi- 
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channel LDV interface module to synchronize the data 
acquisition Q. A triggering pulse from the LDV signal 
processor initiates the acquistion of an analog temper- 
ature signal. The spatial separation between the LDV 
focusing spot and the thermistor tip is kept at a mini- 
mal value of 0.7 ± 0.2 mm, a distance 3 times larger than 
the tip diameter of the thermistor but 20 times smaller 
than the correlation length between the temperature and 
velocity fluctuations 0- The simultaneous velocity and 
temperature measurements can thus be assumed to be 
taken at the same location. In this work, we study mea- 
surements taken at two locations within the convection 
cell: the cell center and near the sidewall (on the mid- 
plane, at 8 mm from the sidewall). Near the sidewall, the 
vertical velocity component v z and one horizontal veloc- 
ity component v y , which is out of the rotation plane of the 
mean large-scale circulation, were measured. At the cell 
center, all three velocity components v x , v y and v z were 
measured with v x along the direction of the large-scale 
mean circulation near the bottom plate. The number of 
velocity or temperature measurements at each location 
is 2 — 4 x 10 5 , corresponding to 8 hours of real-time mea- 
surements. 

From the simultaneous velocity and temperature mea- 
surements, we calculate (v|T(£)), the conditional aver- 
age of the velocity on the temperature measured at time 
t. In calculating (v|T(£)), we take the average of those 
velocity measurements with the corresponding tempera- 
ture measurements falling within T(£) ± 0.005°C. Then 
we decompose the velocity measurement v(£) into a sum 
of (v|T(£)) and the remaining part, denoted as Vb(£): 



v(i) = (v|T(t))+Vb(t) 



(4) 



By construction, Vb(£) is uncorrelated with T(t) and av- 
erages to zero over time. Thus we take it as the back- 
ground velocity fluctuation. If v(£) correlates with T(£), 
(v|T(£)) is different from the ordinary average veloc- 
ity (v) and is a function of T as well as £ since T(£) 
depends on £. In this case, we define the velocity of 
the plume by v p (£) = (v|T(£)). In other words, v p (£) 
is the part of v(£) that correlates with T(£). If v(£) 
is uncorrelated with T(t), Vb(£) = v(£) — (v) and no 
plume is present. With this velocity decomposition, 
we can easily identify the local heat flux carried by 
the plumes. The normalized local heat flux is given 
by j(£) = [v(£)c$T(£)]L/(kA). The part carried by the 



plumes is thus j p (£) 



,(£)5T(£)]L/(kA). When the 



contribution by the plumes is subtracted, the remaining 
is the heat flux carried by the background velocity fluctu- 
ation, which is jb(£) = [vb(£)5T(£)]X/(/«A) and averages 
to zero over time. 

It is found that the horizontal velocity components 
have little correlation with the temperature fluctuation, 
having a correlation coefficient of 0.1 and 0.2 respec- 
tively at the cell center and near the sidewall. On the 
other hand, the vertical velocity component v z shows 



stronger correlation with T(t), having a correlation coef- 
ficient of 0.3 and 0.5 respectively at the cell center and 
near the sidewall. Consistently, we find that the differ- 
ence between (v|T(£)) and (v) lies mainly in the vertical 



direction. That is, v p (£) 



i-'pz 



(£)z. 



It was recently found |7j that the probability distri- 
bution of the local heat flux in the horizontal direction 
is approximately symmetric whereas that in the verti- 
cal direction is skewed towards the positive fluctuation. 
This asymmetry in the vertical flux distribution was in- 
terpreted to be caused by the plumes0. We use these 
results to test our scheme of plume extraction. We expect 
that after subtracting the the contribution by the plumes, 
the probability distributions of the local heat flux in the 
horizontal and vertical directions would become identi- 
cal. In Fig. we show the probability distributions of 
jbxt jby, and jb z . Here the standardized variable X is de- 
fined as (X — (X) ) /<tx where (X) and ax are respectively 
the mean and standard deviation of X. In the calcula- 
tion, we estimate the mean temperature of the bulk fluid 
T by the average temperature measured at the cell cen- 
ter. As expected, the standardized distributions of the 
horizontal and vertical components of jb are the same. 
Moreover, they are approximately the same at the two 
locations studied. 

In a recent studv[loj|. the three velocity components 
at the cell center were found to possess the She-Leveque 
hierarchical structure [Tlj: 



'p+2 
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S P (r) 



(5) 



where S p (t) = (\v(t + r) — v(t)\ p ) with v being one of the 
three velocity components, A p are constants independent 
of r and S^°°'(t) = lim p _>oo S p+ i(t)/S p (t). Here, is a 
parameter whose value lies between and 1. The smaller 
[3 is, the more intermittent the velocity fluctuation is. It 
was found that the two horizontal velocity components 
are characterized by the same value of (3 but the vertical 
velocity component is characterize d by a smaller value of 
(3. This distinction was attributed to the presence of 
the plumes which makes the vertical velocity component 
more intermittent. As a further test of our scheme of 
plume extraction, we analyze Vb(£) at the cell center to 
check whether its horizontal and vertical components will 
now be characterized by the same value of (3. To get the 
value of p, we follow the procedure developed in Ref . 0] 
and calculate the relative scaling exponents p{p) of the 
normalized structure functions, defined by 



S P (r) 
S 2 (t)p/ 2 



Si(r) 



pip) 



(6) 



If Eq. © holds, we have G3 Ap(p + 1) = f3Ap(p) - 
(1 + (3), where Ap(p) = p(p + 1) — p(p). In Fig. [21 we 
plot Ap(p + 1) vs Ap(p). It can be seen that the data 
points for the three velocity components of Vb (£) nearly 
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FIG. 1: The standardized probability density functions 
P(jbx) (triangles), P(jby) (squares), and P(jbz) (circles) at 
Ra = 2.6 x 10 9 . The measurements were made at the cell 
center (top panel) and near the sidewall (bottom panel). 
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FIG. 2: Ap(p+1) vs Ap(p) for the :r-component (triangles), y- 
component (squares), and ^-component (circles) for v (open 
symbols) and Vb (filled symbols) at the cell center for Ra 
= 4.8 x 10 9 . 
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fall on the same line, showing that the difference in the (3- 
values of the horizontal and vertical components indeed 
becomes vanishingly small. 

Next we use our scheme of plume extraction to study 
the temperature dependence of the plume velocity. The 
results are shown in Fig. [21 We find that v pz = aST at 
the cell center and v pz = by/T - T s for T > T (> T s ) 
near the sidewall for some constants a and b. The fitted 
value of T s differs from the value of To as estimated by 
the average temperature at the cell center by about 1%. 
We believe that this difference is caused by drifts in the 
mean temperature of the bulk fluid over time and that 
the measurements at the two locations were not taken at 
the same time. Hence, it is reasonable to conclude that 
v pz oc V ST for ST > near the sidewall. 

In the following, we shall understand this temperature 
dependence. At the cell center, the mean velocity van- 
ishes and thus we expect the viscous term to dominate 



FIG. 3: v pz vs T at the cell center (circles) and near the 
sidewall (squares) at Ra = 2.6 x 10 9 . The solid lines from 
bottom to top are respectively the fits of aST and by/T — T s 
for some constants a and b. The fitted value for T s is 29.9° C 
while T is 30.2°C. 

over the advection term. By balancing the viscous term 
by the buoyancy term in Eq. we have 

center : = gaST (7) 

l c 

for some length scale l c , which explains the temperature 
dependence observed at the cell center. Since plumes are 
excitations from the thermal boundary layer, we expect 
l c to be of the order of the thermal boundary layer thick- 
ness X t h- We calculate l c from the fitted value of a and 
compare l c /L with At^/L^^. As shown in Fig.QJ l c /L, 
after scaled down by a factor of about 2, coincides with 
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FIG. 4: A th /L (triangles) [H and l c /(cL) (filled circles) vs 
Ra with c = 1.9. 



Xth/L over the limited Ra range studied. 

Near the sidewall, the average velocity of the plumes is 
equal to the average vertical velocity, which is essentially 
the mean large-scale circulation velocity and is relatively 
large. Thus the advection term should dominate over the 
viscous term. For ST > 0, we equate the advection term 
by the buoyancy term in Eq. Q and get 



In summary, we have presented a scheme to extract in- 
formation about plumes, the prominent coherent struc- 
tures in turbulent thermal convection. Our scheme in- 
volves a decomposition of the local velocity into two 
parts. The part that correlates with the local temper- 
ature measured at the same time is defined as the plume 
velocity. Using this scheme of plume extraction, we have 
found the temperature dependence of the plume velocity 
at the center and near the sidewall of the convection cell, 
and understood such dependence from the equations of 
motion. Moreover, we have obtained the average local 
heat flux in the vertical direction at the cell center, and 
found that it has a scaling dependence with Ra different 
from that of the Nusselt number. This difference shows 
that heat is not mainly transported through the central 
region but instead through regions near the sidewalls of 
the convection cell. Further results about the mean large- 
scale circulation velocity and the velocity fluctuation at 
the cell center can be derived from Eqs. (0 and (JHJ, and 
will be reported elsewhere. 
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sidewall : 



-f^=gaST 



for ST > 



(8) 



for some other length scale l s . Hence v pz oc V ST for ST > 
reproducing the temperature dependence observed near 
the sidewall. Since the advection term is produced by 
the mean large-scale circulation, it seems reasonable to 
expect l s to be of the order of the cell height L. Indeed l s , 
calculated from the fitted value of b, is about 0.3 — 0.4L 
and has little dependence on Ra. 

We further obtain the average local heat flux in the 
vertical direction at the cell center, (j z (t)) c , which can 
be thought of as a pointwise Nusselt number at the 
cell center. Since (jb(*)) = 0, (j z {t)) = (J pz (t)) = 
(v pz (t)ST(t)}L/(K,A). Using 0, we have 



&(*)>c 



Ra 



((ST) 2 ) c (I 
A 2 



(9) 



Note that A c = yj ({8T) 2 ) C is the root-mean-square tem- 
perature fluctuation at the cell center. It was found 
that A c /A ~ Ra -1//7 both in low temperature helium 
gasQ and in water (lj). As LJL goes like Xth/L (see 
Fig. and X th /L ~ Ra" 2 / 7 0, we have (j z {t)) c ~ 
Ra 1 / 7 . The scaling exponent 1/7 is clearly different from 
the scaling exponent of approximately 2/7 of the Nusselt 
number [lq. This result thus shows explicitly that heat 
cannot be mainly transported through the central region 
of the convection cell. 
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